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Abstract 



pendent operators. In two cases one can construct a state-dependent 
Hamiltonian which is such that the corresponding non-linear von Neu- 
mann equation is solved by the given density operator. We prove that 
in a certain context these two cases are the only possibilities to obtain 
special solutions of this kind. Well-known solutions of the reduced 
Maxwell-Bloch equations produce examples of each of the two cases. 
Also known solutions of the non-linear von Neumann equation in di- 



5_i ■ mension 3 are reproduced by the present approach. 

1 Introduction 

Exact solutions of non-linear differential equations are often based on sets 
of special functions like tanh and sech or Jacobi's double periodic functions 
sn, cn, and dn. Here, the algebraic relations between these functions are 
exploited to construct solutions of non-linear von Neumann equations. Given 
a Hamilton operator H, the von Neumann equation reads 

ip= H,p] . (1) 



1 



The notations 



d 



and 



H,p 



Hp- pH 



(2) 



P 




are used. A non-linear von Neumann equation is obtained when H is allowed 
to depend on p 



The notation H[p] is used instead of H(p) to stress that H[p] is not a function 
of p in the sense of spectral theory. Alternatively, a function f(x) is given 
and the equation is written as ip = [H,f(p)}_. The prototype of the latter, 
obtained with the choice f(x) = x 2 , is the Euler top equation flj and can be 
written as well in the form (JHJ) 



Some explicit solutions [HEj of this equation depend on time t as sech(a;t). 
They have been called [2] self-scattering solutions because they exhibit dif- 
ferent behaviour in the t — > — oo and t — ► +00 limits and make the transition 
from one behaviour to the other in a limited region of time. However, up to 
now, the reason why the sech function appears had not been investigated in 
a systematic manner. 

It has been observed |3j that solutions of the von Neumann equation with 
time-dependent Hamiltonian are also solutions of a non-linear von Neumann 
equation. The main application of the present work supports this point of 
view. The reduced Maxwell-Bloch equations describe the time evolution of 
a spin variable in presence of an applied field. It was noted i4J that the 
well-known McCall-Hahn Ej solution of these equations is a solution of a 
non-linear von Neumann equation as well. This link between time-dependent 
and state-dependent Hamiltonians is clarified in the present paper. 

The non-linear von Neumann equation has mainly been investigated with 
Darboux transformations, trying to find non-trivial solutions (see e.g. 121 
ElEl El El CD!) Many of the solutions constructed in this way concern 3-by-3 
matrices. The connection between these special solutions and Jacobi's double 
periodic functions was noticed in the Appendix of |2j. Here, we reproduce 
these solutions as special cases of our approach. 

In the next section, some properties of solutions of non-linear von Neu- 
mann equations are discussed. In Section 3 we review basic properties of 



ip = H[p},p 



(3) 



ip = H, p 2 



{H,p},p 



with {H, p} = Hp + pH. 



(4) 
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Jacobi's double periodic functions and formulate methods to construct state- 
dependent Hamiltonians for which the corresponding non-linear von Neu- 
mann equations have special solutions. Two cases are formulated and for 
each theorem is formulated which shows that this case is the only 

possible solution of the non-linear von Neumann equation that has the given 
form. Section 4 treats different applications. The paper ends with a short 
summary in Section 5. The proofs of the two theorems are given in Appendix. 



2 Non-linear von Neumann equation 
2.1 Definition 

We restrict ourselves to non-linear von Neumann equations of the form (J3J), 
with a Hamiltonian which depends linearly on the density operator p. E.g., 
H [p] might be of the form 

H[p] = J2^kX*pX k . (5) 

jk 

The operators Xj and the matrix coefficients Xjk do not depend on time. 



2.2 Spectrum 

A nice property of solutions of the non-linear von Neumann equation is con- 
servation of spectrum. Indeed, let be given a solution p t of (JHJ), and let U t be 
solution of 

^U t = iU t H\pt] (6) 
with initial condition Uq = I. Then one has 

Pt = u; Po u t . (7) 

The adjoint operators U£ are isometric. To see this, note that 

^ UtU* t = (8) 

because H[p t ] is self-adjoint. Hence one has UtU* = I, i.e. U£ is isometric. 
Assume now poip = Then there follows ptU*ip = XU^ip. Conversely, 
if p t ij) = Xij) then poU t ip = \U t ip follows. Hence, p t and p have the same 
spectrum. 

Note that the solution of (jSJ) can be written as 

U t = I + i dt 1 #[p tl ]+* 2 / dti / dt 2 H[p t2 ]H[ Ptl } + --- . (9) 
Jo Jo Jo 
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2.3 Linear part 

The Hamiltonian H[p] may contain a contribution of the form (Tr p)Hq, 
which for properly normalized density operators leads to the standard linear 
von Neumann equations. For simplicity, this contribution is called here the 
linear part. Quite often it can be eliminated, leading to a problem which can 
be discussed more easily. 

Let p t be a density operator satisfying the non-linear von Neumann equa- 
tion (JHJ). Let Hq be a symmetric operator satisfying for all t 

H[e itHo p t e- itHo ] =e itHo H[p t ]e- itHo . (10) 

Introduce a linear map K[a] by 

K[a] = H[a]-(Tra)H Q . (11) 

Then a t , defined by 

a t = e UHo p t e- UHo , (12) 
satifies the non-linear von Neumann equation (JHJ) with H [p] replaced by K[a] 

ta t = [K[a t ],a t ]_. (13) 



3 Special solutions 

3.1 Jacobi's elliptic functions 

Solutions of non-linear equations are often based on Jacobi's elliptic functions 
with elliptic modulus k. They satisfy the algebraic relations 

1 = sn 2 (x, k) + cn 2 (x, k) 
1 = dn 2 (x, k) + k 2 sn 2 (fc, x) 

— sn(x, k) = cn(x, k) dn(x, k) 
d 

— cn(x, k) = — sn(a;, k) dn(x, k) 
d 

— dn(x, k) = —k 2 sn(x,k)cia(x,k). (14) 
ax 

In the limit k = 1 the period of these elliptic functions diverges and sn(x, k) 
converges to tanh(x), cn(x, k) and dn(x,k) both converge to sech(x). The 
above relations then simplify to 

1 = sech 2 u + tanh 2 u 
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— sechu = — sechw tanh u 
d 

— tanhw = sech 2 w. (15) 
du 

In the limit k = the function sn(x,k) converges to sin(x), cn(x,k) to 
cos(x), and dn(x, k) converges to 1. 

We next consider special solutions of the non-linear von Neumann equa- 
tion involving three independent operators. Two different cases are consid- 
ered. 



3.2 Case 1 

Theorem 1 Let be given p(t) of the form 

pit) = 8+ cn(ujt, k)A + sn(ujt, k)B + dn(ut, k)X. 
Assume pit) satisfies the nonlinear von Neumann equation 

P = 1 [p,H[p]]_ . 

with H[-} linear. 

Assume that A, B, and X linearly independent operators. 
Assume that 6 commutes with A, B , and X . 
Assume uj ^ and < k < 1. 

Assume that the range of H[-] is in the linear span of the operators A, B, 
and X . 

Then there exist constants a and (3 such that 

i[B,X] = aA (16) 

i[A,B] = k 2 (3X (17) 
ad 

t[A,X] = — Q B. (18) 



a + /3 



The Hamiltonian is of the form 



H[A] = (v+^y (19) 

H[B] = uB (20) 

H[X] = (u-^)X (21) 

H[9] = 0. (22) 

Conversely, given operators A, B, X, and 9, satisfying the above assump- 
tions, then the operator p(t) given by Mb}) solves the non-linear von Neumann 
equation. 
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One direction of the proof is given in Appendix A. The proof of the 
converse statement is a straightforward exercise. 



3.3 Case 2 

Theorem 2 Let be given p(t) of the form 

p(t) = 9 + cn(ut, k)A + sn(ut, k) da(ut, k)C + cn(ut, k) 2 D, 

with 9, A, C , and D linearly independent operators, and with uj ^ and 
< k < 1. 

Assume p(t) satisfies the nonlinear von Neumann equation 

p = i [p,H[p]]_ . 

with if [•] linear. Assume that the range of H[-] is in the linear span of the 
operators A, C , and D. 

Assume that an operator # exists, which commutes with A, C , and D, 
such that 9 — 6 belongs to the linear span of A, C, and D. 

Then there exist constants a and 5 such that 

t[C,D] = aA (23) 
i[A,C] = 5D (24) 
i[A,D] = -k 2 5C. (25) 

The Hamiltonian is necessarily of the form 

H[A] = (u + ^)a (26) 



S f 

H[C] = uC (27) 
H[D] = uD (28) 

H m - {~ + 1 -^f" a+ » D - < 29 > 

The operator 8q is given by 

Conversely, given operators A, C , D, 9 and 9, satisfying the above assump- 
tions, then the operator p(t) given by Q2fy) solves the non-linear von Neumann 
equation. 

The difficult direction of the proof is given in Appendix B. The proof of 
the converse statement is a straightforward exercise. 
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4 Examples 



4.1 Reduced Maxwell-Bloch equations 

Let oi, 02, 03 denote the Pauli matrices. Consider the density operator 



3 



1/-, . x 1 / 1 + % Ui-%u 2 \ , Q1 x 

p t = -(1 + > u a a a ) = - (31) 

a—1 



with 



2 

M 2 = ; : f 7 x 9 sech(t/r)tanh(t/r) 
1 + (rAj^ 

"a = -1 + 1 - ( 2 rA)2 sech 2 (t/r) . (32) 

See (Eq. (4.21) of 0). It solves the reduced Maxwell-Bloch equations 

111 = —Au2 

112 = Ami + k£v,3 

U3 = —k£u2 (33) 
with constants A, k, and r and with 

Sit) = — sech(t/r). (34) 

KT 

This solution satisfies the requirement that \u\ is constant in time so that 
the eigenvalues of p, which equal (1/2)(1 ± \u\), are conserved. In fact, one 
has \u\ = 1. Hence, p is a one-dimensional projection operator. 

This density operator is an example of Case 2, with k — 1, to — 1/r, and 
with operators 

= \l (35) 



rA 
1+~^A) J 



A = i I (36) 



c ^ TTW ff2 (37) 
D = TTW^ (38) 
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They satisfy the commutation relations 

2rA 

' lAC] - = " TTW C (39) 

! i c - c i- = -^rrw- 4 (40) 

2rA 

!|fl " 41 - = -TTW C ' (41) 

From Theorem 2 follows that pt satisfies the non-linear von Neumann 
equation with Hamiltonian given by 

H[a- L ] = -I^ + AV (42) 

H[a 2 ] = (43) 

H[a 3 ] = (44) 

9/ .3 

ff M = "Z?+A^- (45) 

4.2 Including phase modulation 

A related problem is obtained by including phase modulation effects. The 
equations are 

ill = 4>U2 

U2 = — <Mi + kSus 

u 3 = —n£ii2 (46) 



with 



KT 

Solutions are (Eq. (4.49) of 0) 

U\ = t5u2 



(p = — 5tanh(£/r) 

£ = + (r5) 2 sech(t/r). (47) 



U2 = -7rTw sech(t/r) 

u 3 = tanh(t/r). (4E 



Note that 



1 

U\ = M1M3 

T 
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1 

u 2 = — u 2 u 3 

T 

us = -(l-ul) = \u\ + ul). (49) 

r r 

Let us try to match this example to Case 1. Clearly needed is u — 1/r, 
k 2 = 1, 9 = (1/2)1, and B = (1/2)03. The choice of operators A and X is 
not so obvious. But the choice 

1 t5 

A = -- - = = 0x (50) 

2 VI + rW V ^ 

X = -- - 1 =0 2 (51) 
2 VI + r 2 5 2 1 ^ 

(52) 

satisfies the requirements. From Theorem 1 then follows that p t satisfies the 
non-linear von Neumann equation with Hamiltonian given by 

#[<n] = (53) 


H[a 2 ] = -5a 2 (54) 

H[a 3 ] = (55) 

H[I] = 0. (56) 

4.3 Three-level system 

Fix < k < 1 and a, 5, and <fi real. Consider 3-by-3 matrices 

e i4, 

A = k^= I e*+ I , (57) 

-ie l 

C = x /^- ( ie* I , (58) 

D = kVa5 [0-10], (59) 
\ / 

0o = h- (60) 
















e -i< 
















ie~ 


i<t> _i e - 




' 





-1 









It is easy to verify that they satisfy the case 2 commutation relations 

. Hence the density matrix p(t), given by (|23|). satisfies the non-linear von 
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Neumann equation (fT6*j) for any Hamiltonian of the form (|29|) . The operator 
9 is given by (f3*Uj). 

This result can be transformed by adding a linear part (in the sense of 
Section 2.3). Let a{t) = e'^ 3 p{t)e ifitPi with 



P 3 = . (61) 



Note that 

H[a(t)} = e-^ tPz H{p(t)}e^ m (62) 

is trivially satisfied because H[A] is proportional to A, H[C] is proportional 
to C, H[D] is proportional to D, and 9 commutes with P 3 . This implies that 
a(t) satisfies the equation 

id = [H[a] + /2P 3 , a] (63) 
The latter equation can be written into the form 

id = [H + Hj, a] , (64) 

where 



i*o = -A , (65) 



and 






Hi = e m(ut, k) [ e^~^ | , (66) 

e ~i(<j>-fit) e -i{4>-fjLt) 




and A = —ku^J^, e = v = 0. It describes a three-level system interacting 
with an electromagnetic pulse E = Eoe % ^~^ cu(ut, fc). 
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4.4 Known solutions in d = 3 

Fix real constants k, u, 0, A and /i, satisfying |A| < /i and < k < 1. Define 
operators 9, A, B, and X, by 

= \l, (67) 
^ = I 1 ] (68) 




A = (69) 



X = -ie** . (70) 



They satisfy the case-1 commutation relations with 

" = ^ < 72 > 
Hence, the results of the previous section imply that the density operator 

a t = 9 + B sn(wt, fc) + A cn(ut, k) + X dn(wt, fc) (73) 

satisfies the non-linear von Neumann equation «cr = cr]_ for any Hamil- 

tonian -f/"[cr] satisfying 

(74) 
(75) 
(76) 
(77) 

A satisfactory choice is 



(78) 



H[9] = 





H[A\ = 


M 


H[B\ = 




H[X] = 







0"31 





-0-32 


Cl3 ~ °23 












0"31 








0"32 


0"13 


0"23 






Note that 



H[a] = {H , a} -^(Tr a)H (79) 
11 




with 

H = ( -fi | . (80) 
In particular, p t defined by 

Pt = e -(V^tH 0(Tte (2/3)UH (gl) 

satisfies the non-linear von Neumann equation ip = [Ho,p 2 ]_. This equa- 
tion with Hamiltonian (jSOJ) has been studied in [2|. Explicit solutions were 
obtained in the limit k = 1. 



4.5 Variations on a theme 

Let us slightly modify the previous example. Fix real constants b ^ 0, u, <f>. 
and < k < 1. Define operators 9, A, B, and X, by 

9 = -I, (82) 



A = —= 1 (83) 




B = — (84) 




X = -^f= -ie^ I . (85) 



bV2 



ie-V 



They satisfy the case-1 commutation relations with 

a = /? = --. (86) 

The density operator 

a t = 6 + B sn(wt, fc) + A cn(wt, fc) + X dn(wt, fc) (87) 

satisfies the non-linear von Neumann equation zct = <r]_ for any Hamil- 
tonian i?[cr] satisfying 

ff[0] = (88) 

H[A] = (u-b)A (89) 

H[B] = vB (90) 

H[X] = (u + b)X. (91) 
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A satisfactory choice, corresponding with v = 46, is 



With this choice is 



with 



3(T 2 1 4(731 

H[a] = b ( 3(7x3 5(T 32 J . (92) 

4(7 13 5(723 



if[<r] = {i/o^}-^(Tr(T)//o, (93) 



1 

H = b ( 2 | . (94) 
x 3 

In particular, p 4 defined by 

pt = e -(2/3)itH 0(Tte {2/3)itH (Q5) 

satisfies the non-linear von Neumann equation ip = [Ho,p 2 ]^. In the limit 
k — 1 this example has been discussed in [I]. 

5 Discussion 

This paper studies the non-linear von Neumann equation under the restric- 
tions that (1) the Hamiltonian H[p] depends linearly on the density operator 
p; (2) the solution p t involves Jacobi's elliptic functions sn, cn and dn, or 
there limits tanh and sech; (3) the solution p t involves 3 linearly indepen- 
dent operators, and possibly a fourth operator commuting with these three 
operators. 

The paper does not focus on methods for solving non-linear von Neu- 
mann equations. It rather follows the opposite way. Starting from a special 
solution p t , it constructs the Hamiltonian H[p] for which p t solves the corre- 
sponding non-linear von Neumann equation. Two cases have been treated. 
In both cases we have proved a theorem stating conditions under which the 
special solution is found. Two well-known solutions of the reduced Maxwell- 
Bloch equations are examples of these two cases. We have shown that the 
reduced Maxwell-Bloch equations can be generalized to three-level systems 
and that these have not only special solutions involving sech and tanh, but 
also periodic solutions involving Jacobi's elliptic functions. Finally, some of 
the known 3-dimensional solutions of the non-linear von Neumann equation 
appear to be examples of these two cases as well. 
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The examples of the reduced Maxwell-Bloch equations show that non- 
linear von Neumann equations and their solutions appear in physics in a 
natural manner. The d = 3-solutions of the equation ip = [Ho,p 2 ], obtained 
here, are generalizations of those found in 0j. 
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A Proof of Theorem 1 

A straightforward calculation leads to the set of equations 



= i[H{A},A]_+i{H{X},X]_ (96) 

= -i[H[A],A]„ + i[H[B],B}_ - k 2 i[H[X],X}_ (97) 

= i[#[0],A]_ (98) 

= i[H[6],B]- (99) 

= i[H[9],X\- (100) 

uA = i[H[B],X]-+i[H[X\,B]- (101) 

-coB = i[H[A],X]-+i[H[X],A]- (102) 

k 2 ujX = i[H[A],B]-+i[H[B],A]-. (103) 

Equations JHEJ EH HHOJ) imply that H[9] commutes with A, B, and X. 
Introduce the notations 

H[A] = a A A + b A B + x A X (104) 

H[B] = a B A + b B B + x B X (105) 

H[X] = a x A + b x B + x x X (106) 

H[9] = a A + b B + x X. (107) 
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Then the remaining equations become 



= -b A K AB + (a x -x A )K AX + b x K BX (108) 

= (a B + b A )K AB + {x A - k 2 a x )K AX - (x B + k 2 b x )K BX (109) 

= -b K AB -x K BX (110) 

= a K AB -x K BX (111) 

= a K AX + b K BX (112) 

uA = a x K AB + a B K AX + (b B - x x )K BX (113) 

-ujB = —b x K AB + (a A — x x )K AX + b A K BX (114) 

k 2 uX = (a A - b B )K AB - x B K AX - x A K BX . (115) 

Because A, B, and X are linearly independent and uo and k do not vanish 
the last three equations imply that K AB , K AX , and K BX are linearly inde- 
pendent. Then the first five equations imply a B = a x = b A = b x = x A = 
x B = a = b = c = 0. The three last equations then read 

uA = (b B -x x )K BX (116) 
-loB = (a A -x x )K AX (117) 
k 2 uX = (a A -b B )K AB . (118) 



This implies (JH EH HHJ) with 



a = (119) 



Let v = b B . Then one finds 



The equations for H[-} then follow. 



b B - x x 
aA — b B 



P = -1T-T-- ( 12 °) 



a A = (121) 
x x = u-~. (122) 



a 



15 



B Proof of Theorem 2 



B.l Equations 

A straightforward calculation leads to the set of equations 






= i[D,H[D]}-k 2 i[C,H[C}} 


(123) 





= i[9,H[9]] + (l-k 2 )i[C,H[C]] 


(124) 





= i[9,H[D]}+i[D,H[9))+i[A,H[A))- 


f (2k 2 -l)i[C,H{C\\2h) 





= i[C,H[D])+i[D,H[C]) 


(126) 


-uoA 


= i[O,H[C\]+i[C,H[0\] 


(127) 


2k 2 uC 


= i[A,H[D]]+i[D,H[A]) 


(128) 


uC 


= i[9,H[A)]+i[A,H[9)}+i[A,H[D}}^ 


-i[D,H[A\] (129) 


-2coD 


= i[A,H[C]]+i[C,H[A}}. 


(130) 



Introduce the notations K A b = i[A,B]-, and similar notation for other 
commutators. Because the range of H[-] is in the span of A, C, and D, one 
can write 



H[A] = a A A + c A C + d A D (131) 

H[C] = a c A + c c C + d c D (132) 

H[D) = a D A + c D C + d D D (133) 

H[0] = a A + c C + d D. (134) 

By assumption there exist numbers t A , tc, and t^, and an operator 9 com- 
muting with A, C, and D, such that 

9 = 9 + t A A + t c C + t D D. (135) 

This implies 

K M = t c K AC + t D K AD (136) 

K ce = -t A K AC + t D K CD (137) 

K m = -t A K AD -t c K CD . (138) 

Then the equations ()127| I128pi3()|) become 

- uA = (c c t A - a c t c - a )K AC + (d c t A - a c t D )K AD 

+(d c t c - c c t D + d )K CD (139) 

2k 2 uC = c D K AC + (d D - a A )K AD - c A K CD (140) 

-2ujD = (cc - a A )K AC + d c K AD + d A K CD . (141) 
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Because A, C and D, are linearly independent and k and 00 do not vanish one 
concludes that the operators Kac-i Kad, and K C d, are linearly independent. 
Consider next (jl23[ I126J) . They can be written as 



= -(c D + k 2 d c )K CD - a D K AD + k 2 a c K AC (142) 
= (d D - c c )K CD - a c K AD - a D K Ac . (143) 

Because of the linear independence of K AC , K AD , and K C d, they imply 
etc = cld = and Cd = —k 2 dc and do = cc- The equations (|139f 114()| 1141ft 
can therefore be written as 



-uA = {pet a - ao)K AC + d c t A K AD 

+(d c t c - c c t D + d )K CD (144) 

2k 2 uC = -k 2 d c K AC + (c c - a A )K AD - c A K CD (145) 

-2uD = {c c - a A )K AC + d c K AD + d A K CD . (146) 

Finally, (fHHjl . (fT^Kjl. and (TT^l. the latter simplified with (fT^Hjl . become 

= (c tA - a t c )K AC + (d t A ~ a t D )K A D 

+ (d t c - c t D + (1 - k 2 )d c )K CD (147) 

= (ca - k 2 d c t A )K A c + (cctA + d A - ciq)K A d 

+ (k 2 d c t D + c c t c + (2k 2 - l)d c - c )K CD (148) 

= (-k 2 (l-2k 2 )dc-2k 2 (c -a A tc + CAt A ))K A c 

+ ((1 - 2k 2 )(c c - a A ) - 2k 2 (d - a A t D + d A t A ))K A D 

+(_(! _ 2k 2 )c A + 2k 2 (c A t D - d A t c ))K CD . (149) 

Because of the independence of Kac, Kad, and Kcd, nine equations follow 

= c t A - a t c (150) 

= d tA - a t D (151) 

= d t c - c t D + (1 - k 2 )d c (152) 

= CA -k 2 d c t A (153) 

= c c tA + d A -a (154) 

= k 2 d c t D + c c t c + (2k 2 - l)d c - c (155) 

= -k 2 (l - 2k 2 )dc-2k 2 (co-a A tc + c A t A ) (156) 

= (l-2A; 2 )( Cc; -aA)-2A; 2 (rfo-aAtD + ^) (157) 

= -(l-2k 2 )c A + 2k 2 (c A tD- d A t c )- (158) 



These equations fix a , Co, do, ca, and in terms of the remaining parame- 
ters. The latter are constraint by another 4 equations. Their analysis is done 
in the next subsection. 
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B.2 Analysis 

Lemma 1 ca = 0. 

Proof Multiply fT^j) to obtain 

(1 — k 2 )a dc = doC t_D — aodot c . (159) 

Then (115011151]) imply that the r.h.s. vanishes. One concludes that a d c = 0. 
This leaves two cases 

1) dc = 0; the lemma follows from (|153j) . 

2) dc 7^ and ao = 0; then ()152)) implies that cq and do cannot both vanish. 

But this implies = via (jl5(ty I151j) . Again ca = follows from 




Lemma 2 Either one of the following cases holds 

• a) dA = and a^dc = dctA = 0; 

• b) dA 7^ and tc = dc = Co = 0. 

Proof From ()158j) . using = 0, follows d^tc = 0. Hence either dA = or 
dA 7^ and = 0. 

Use that a^dc = dctA = from the proof of the previous lemma. Then 
()154j) . multiplied with dc, implies that dAdc = 0. Hence the assumption 
that dA 7^ implies dc = 0. But 4^0 implies also tc = 0. Finally, ()155|) 
implies c = 0. 

Next assume 0^ = 0. The results a d c = dctA = are taken from the 
proof of the previous lemma. 

□ 

The first possibility of the previous lemma is the only one that can occur. 
This is proved in following lemma. 

Lemma 3 dA = 0. 

Proof The proof goes ex absurdo. Assume that the second case of the 
previous lemma holds, i.e. dA ^ and ca = tc = dc = c = 0. Then the 
operator expressions ()144| I145| 1146}) become 

-uA = (c c tA-ao)KAc + (-cctD + do)K C D (160) 
2k 2 uC = (cc-a A )K AD (161) 
-2uoD = (c c - a A )K AC + d A K CD . (162) 
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The solution of this set of equations is 

K AC = -[-d A uA + 2(d -c c t D )uD) (163) 
2k 2 

K AD = uC (164) 

c c - a A 

K C d = [{cc ~ a A )uA - 2(c c t A - a )uD] , (165) 

with 

N = d A (c c t A - do) - {cc ~ a A )(d - c c t D ). (166) 
The Jacobi identity then implies 

= i[A,Kc D )-i[C,K AD )+i[D,K A c} 

= jj[-2(c c t A -a )+d A )uK AD . (167) 

Hence there follows d A = 2(cct A — do). Comparison with ()154j) then implies 
d A = 0, which contradicts the assumption that d A ^ 0. 

□ 

Lemma 4 d c = 0. 

Proof From Lemma El follows dct A = 0. This implies dc = or t A = 0. Let 
us assume t A = 0. Then, using the results of the previous lemmas, Q150I - 
I158j) reduce to 

= d t c - c t D + (1 - k 2 )d c (168) 
= a Q (169) 
= k 2 d c t D + c c t c + (2k 2 - l)d c - c (170) 
= (l-2k 2 )d c + 2(c -a A t c ) (171) 
= (l-2k 2 )(c c -a A )-2k 2 (d -a A t D ). (172) 

Use the last two equations to eliminate Co and do from the first and the third 
equation. This gives 

= (1 - 2k 2 ) [k 2 d c t D + (c c - a A )t c \ + 2k 2 (l - k 2 )d c (173) 
= k 2 d c t D + (c c - a A )t c + \(2k 2 - l)d c . (174) 
Combining these two equations gives dc = 0. 

□ 
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Lemma 5 c = t c = 0. 



Proof From (j!45J) . using ca = d c = 0, follows cq ^ cla- But from ()1551 1156]) 
follows 

= (c c - a A )t c . (175) 
Hence, tc = follows. From ()155|) or ()156j) then follows Co = 0. 

□ 



Lemma 6 = t A = 0. 

Proof From (TED EH follows 

= (do - cc*^. (176) 

Note that <io 7^ c^to follows from ()144j) . One concludes therefore that t A = 0. 

Finally, ^4 = and ()154|) imply cc-t^ — a = 0. Hence t A = implies 
ao = 0. 

□ 



B.3 Completing the proof 

Thus far, the desired form of the operators in (J2~31 — IHUj) has been proved. The 
actual relations between the coefficients follow by inserting these expressions 
into the equations (jl2ril — [TT?UJ) . In particular, the operator expressions (JHU 
11351 Eg]) become 

2fc 2 ^C = (c c -a A )K AD (177) 
-2^D = (c c -a A )K AC (178) 
a;A = (c c t D - d )K CD . (179) 



The solution of this set of equations is of the form (|2*Hl l2~4l l23j) with 

2cu 



Cc - CLA 



;iso) 



a = ^— . (18r 
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